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The mean time for a diffusing ligand to activate a target protein located on the surface of a 
microdomain can regulate cellular signaling. When the ligand switches between various states 
induced by chemical interactions or conformational changes, while target activation occurs in only 
one state, this activation time is affected. We investigate this dynamics using new equations for the 
sojourn times spent in each state. For two states, we obtain exact solutions in dimension one, and 
asymptotic ones confirmed by Brownian simulations in dimension 3. We find that the activation 
time is quite sensitive to changes of the switching rates, which can be used to modulate signaling. 
Interestingly, our analysis reveals that activation can be fast although the ligand spends most of the 
time 'hidden' in the non-activating state. Finally, we obtain a new formula for the narrow escape 
time in the presence of switching. 



Cellular chemical reactions depend on the activation of small targets by diffusing ligands. This process can be 
regulated by several parameters such as the geometry of the cellular microdomain, the target shape, or the state of 
the target and ligand upon encounter. In the past, the activation rate and survival probability of the target were 
studied for diffusion-influenced chemical reactions where the reactivities of the target or the ligands stochastically 
fluctuate in time [U, Interestingly, for slow gating dynamics it was found that both processes are not equivalent, 

and, for example, a gated target may lead to non exponential behaviour at short times P, When the target 

has a fluctuating potential barrier, there is an optimal combination of parameters for which the mean activation time 
is minimal, leading to a resonant-like phenomena 0, cellular microdomains, target activation may depend on 

the state of the ligand. This is for example the case in the cytoplasm, where enzymes can switch between an inactive 
and active state, or in the nucleus, where a transcription factor needs to be first activated in order to bind to a 
specific DNA promoter ■ Recently, intermittent search scenarios were introduced and extensively analyzed for a 

dynamics switching stochasticallybetween fast ballistic phases and slow diffusive phases, while the target can only be 
found in the diffusing phase [ill. Il2l. [isj: interestingly, switching can decrease the search time, and there are optimal 
search strategies that minimizes the search time. 

In the absence of switching, target activation is determined by the narrow escape time (NET), which is the mean 
time for a Brownian ligand to find a small target in a confined environment [H, [H [H, [13, . The NET has been 
used to compute the mean and variance of chemical reactions with few molecules diffusing in microdomains [l9| . to 
estimate the probability and the arrival time of viral particles to nuclear pores [20| . or to study the early steps of 
phototransduction in rod photoreceptor [2l| . 

We study in this letter the gated narrow escape time (GNET) to exit the domain fl through a small window, if a 
diffusing ligand stochastically switches between two states 1 and 2 with diffusion coefficients Di and D2, and can exit 
only in state 1. Switching may be due to conformational changes or chemical interactions. We estimate the GNET 
using new equations for the sojourn times the ligand spends in the different states. We find that switching not only 
affects drastically the exit time and the sojourn times, but also, only for D2 > Di the GNET can be optimized as a 
function of the switching rates. In addition, a ligand may exit almost as fast as possible, although it spends most of 
the time in state 2. Finally, we give a new formula for the GNET in dimension three, which extends the NET formula 
to the switching case. We also discuss briefiy possible applications in cellular signaling. 

A Brownian ligand diffuses in a confined domain Q. while switching between states 1 and 2 with Poissonian rate 
constants ki2 and k2i and diffusion constants Di and D2. Upon encounter, the ligand activates the target located on 
the small boundary portion dQ,a only in state 1 , while it is reflected otherwise everywhere on the boundary. To study 
the GNET, we consider the mean sojourn times u„(a;,m) the ligand spends in state n before exiting, conditioned 
on starting at position x in state m. From the backward Chapman-Kolmogorov equation [2^ . [23j we find that the 
Un{x, m) satisfy the coupled system of equations 
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L*^Un{x, "l) - ^ krm{Un{x, m) - Un{x, i)) = Snm, (1) 
1=1 

where is the backward Kolmogorov operator in state m, which in our case is L'^ = Dm^, and we have absorbing 
boundary conditions on dVta for u„(a;,l), and reflecting conditions otherwise. Eq. [T] separately constitutes a closed 
system of equations for each value n, and it is sufficient to study the equations for n = 1, because the solutions 
U2{x, ni) are obtained through the linear transformation 

(U2{X, l)\ _ fcl2 A 0^ ^^\^( ^\ (9) 

[u2{x,2) - k2i [O ij Ui(^,2) j + UrJ ■ 



FIG. 1: Example of a trajectory of a diffusing Brownian ligand in a confined domain il and randomly switching between two 
states 1 (continuous line) and 2 (dashed line). In state 2, the ligand is reflected all over the boundary, while in state I it is 
absorbed at dfla- 



The mean sojourn times Ui(l), ui(2), ^2(1) and ^2(2) are obtained by averaging eq. [T] and eq. [5] over an initially 
uniform distribution, satisfying ui(2) = ui(l), ^2(1) — 1*1(1)^12/^1 and ^2(2) — ^2(1) + 1/^21- The mean times 
m(1), u(2) and u to exit the domain starting uniformly distributed in state 1, 2, and in state 1 and 2 with equilibrium 
probability (^1,^2) = (17^, fcllife)' ^re 

U{1) = Ml(l)+U2(l) - Wl(l)(l + fcl2/fc2l) (3) 

u{2) = iti(2) + M2(2) = «(!) + l/fc2i (4) 

u - pMl) + P2u{2) = u{l) + . (5) 

k2i{ki2 + k2i) 

The GNET in one dimension 

We now solve the GNET problem in one dimension where H. reduces to the interval < x < L with an absorbing 
boundary at x = in state 1, and a reflecting boundary at x = L. The onc-dimensional analysis already reveals many 
features that remain valid also in higher dimensions where only asymptotic results are available. Using the scaled 
parameters 

Zi = fci2LV^i, l2 = k2iL'/D2, K = Di/D2 (6) 
we solve eq. [T] and eventually obtain 



ui{x,l) = 



Di {h + l2)Vh+l2 sinhV^i +^2 
I2 (2L — x)x 



h + I2 2Di 
from which we derive for the averaged sojourn time 



(7) 



Ml) = Ti - -A^ (n - ^fih + I2] , (8) 



^1 +^2 V ^1* 

where ri = L'^/{3Di) is the mean first passage time to exit in state 1 without switching, and f{x) = '=°*^^y^ ~ ^ 
monotonically decreasing from 1/3 at cc = towards for x — > 00. In fig. [2 we plot wi(l) as a function of ^i and ^2- 
Interestingly, because wi(l) < ti, we found the non-intuitive result that the sojourn time the ligand spends in state 
1 before exiting is always smaller than the mean first passage time ti to exit in this state without switching. Even 
more surprising, mi(1) can become arbitrarily small by increasing li (resp. the switching rate ^12). This behavior 
can be understood as follows (see also Q): for a ligand that starts uniformly distributed in state 1, the probability 
to find it in the neighborhood of the absorbing boundary at 2: = decreases quickly. But, in state 2 the probability 
distribution is re-homogenized and after switching back to state 1, the density around x = is higher compared to 
the non-switching case leading to a faster exit. Because the sojourn time wi(l) can be arbitrarily small, we wonder 
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FIG. 2: Graph of the sojourn time ?ti(l) obtained from eq. |8]as a function of h and h, scaled by the mean first passage time 
n. In the left panel, the continuous curve is the asymptotic 3/\/li for li ^ 1 and \/Ti 2> h- 



how this affects the mean time to exit when starting initially in state 1. Using eq. [8]we obtain 



"(1) = ^1 + 1^ (^/(^i + '2) + ('^ - l)^^i(l) ) , 



(9) 



which shows that > ti for k > 1. Obviously, by switching to a state with a smaller diffusion constant, we 
cannot speed up exit. However, when k < 1, interestingly, the situation changes: for fixed I2, by expanding 
in eq. [9] as a function of li for li ^ 1, we find that initially decreases if I2 > hi^), where 12{k) is the root of 
f{h) + (k — l)/3 = 0. Together with the asymptotic ~ for large h, we conclude that u{l) has a minimum 
u{l)m < Ti as a function of li for I2 > hif^)- Similar to the results of ^11], we found here that the exit time is 
minimal for a certain value of ki2 (that depends on I2 and k). The left panel of fig[3]shows as a function of li for 
various I2, and the right panel displays as a function of I2 for various h, showing that also has a minimum 
as a function of ^2- For given values I2 and k < 1, we now study the minimum u(l),n of achieved at li ~ h^m- 
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FIG. 3; Graph of the exit time from eq. |9]as a function of h and h for k = 0.1, scaled by n. 

The left panel of fig. [?] shows ui^m scaled by ri as a function of I2 for various k {li.m is found by numerically solving 
du{l)/dli — 0). An asymptotic analysis for large I2 reveals that the position (h^rmh) of the minimum is defined by 



,3/2 

n,m 



3k 



2(1 - k) 



(10) 



from which we derive for large I2 the asymptotic behaviour ui(l)r 
7-1/3 



0(?2 U2{l)r 



T2 



0{q^'^) and u{l)r. 



T2 + 0{l2 ), where T2 = kti is the mean first passage time for a ligand diffusing with a diffusion constant D2. 
The right panel of fig. H] displays the ratio Mi(l)m/'''i and U2(2)„i/t2 as a function of I2 for different k. Surprisingly, 
when exit is fast and u(l) approaches its lower limit T2, the ligand spends almost all its lifetime in state 2 where it 
cannot exit. When considering u(l) as a function of I2 for a given value of li, the minimum position (li, ^2,m) differs 
from {h^rmh), however, it satisfies a very similar relation compared to eq. 1101 In general, for k < 1, u{l) has no 
local minimum as a function of {hyh), but is lower bounded by T2. Only for k > 1 we have the global minimum 
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w(l)m = Tl trivially attained for li = when the ligand does not switch and stays in state 1. Furthermore, for k = 1 
there is a strategy to exit in almost minimal time Ti while spending most of the time hidden in the state 2 where 
exit is not possible. With a 1, the strategy is to choose a^/h ^ 1 and h — ah such that « ri, and since 
Mi(l)/u2(l) — ct, it follows that the ligand spends only the small fraction a/{l + a) of its time in state 1. 

In summary, the lower limit of u{l) corresponds to a ligand diffusing all the time with the maximal diffusion 
constant, and interestingly, fast exit can be achieved even when diffusing most of the time in the state where exit is 
not possible. Furthermore, for k < 1, has a minimum as a function of h resp. I2 for non- vanishing switching rates, 
however, as shown in figure [31 the graph for around and past the minimum is quite flat, while it decays steeply 
for small switching rates. Thus, the behaviour at small rates is an efficient mechanism to modulate the activation 
time, and thus cellular signaling. 




FIG. 4: The left panel displays the minimum u{l)m scaled by ri as a function of I2 and k. The horizontal lines are k, and 
the continuous curves are calculated using eq. 1101 The graphs do not extend up to ^2 = because a minimum exists only 
for I2 > hi^.) (see text). The right panel shows the ratio Ui{l)m/Ti (curves below the horizontal line) and U2(l)m/T2 (curves 
above) of the corresponding sojourn times. The continuous curves are calculated using eq. 1101 

The GNET in a 3- dimensional microdomain 

We now extend our previous analysis to a ligand confined in a 3-dimensional domain U, that activates upon hitting 
a small boundary patch dVla in state 1. Without switching, the time reduces to the narrow escape time (NET) 
0, 0, [i3| ■ The NET analysis showed that outside a small boundary layer around the absorbing hole of radius a few 
g, where a characterizes the extent of cJOa, the positional NET is almost independent of the initial ligand position 
[l3|. Using the scaling introduced in [l^l, we define the dimensionless variables and functions x = x/a,li = ki2a^/Di, 
h = k2ia?/D2, vi{x) = ^^Ui(a7,l), V2{x) = ^^ui(a?,2), and obtain from eq. [T]the scaled equations 

Awi(i;) - li{vi{x) - V2{x)) = 

l^V2{x) + l2{vi{x) - V2{x)) ^ , ^ ' 

where |f2| — \Q,\/a? ^ 1. The boundary conditions are absorbing on dQ,a for wi(a;), and otherwise reflecting. The 
solutions of eq.[Tl]for general values li and I2 are not at hand, thus we present here asymptotical results that clarify the 
effect of switching. For /i ^ 1 or Z2 » at leading orders, vi{x) is solution of the NET problem Awi(a;) — 
When I2 <C ^/h, the leading order solution for vi{x) is found by solving Aui(a;) — l-i{vi{x) — vi) — — where vi 
is the spatial average of vi{x). Using the asymptotic solutions for vi{x), the leading order expression of the sojourn 

time u(l) = (1 + ki2/k2i)^vi is 



Tl , Zi < 1 or /2 > h 

\dna\^/D^2 

To confirm the behavior of u{l) and ui(l) as a function of li and I2, we used Brownian simulations together with the 
Gillespie-algorithm 25] to model switching in a sphere of radius r = 30 with a circular hole of radius a = 1. In the 
left panel of fig.El we show simulation results for mi(1) as a function of li for various Z2. We obtain that ui(l) < ti , 
and confirm the asymptotical behaviour ui(1)/ti « 4/(7r-yi]") which follows from eq. [T^lby using ri w |ri|/(4a_Di) |18| 
and \d^a\ = 7ra^. In the left panel of fig.[5]we display the simulations results for u(l) as a function of h for various I2 
and K — 0.1 {Di — 1, D2 = 10). The plot shows that u(l) has a minimum smaller than ti which is attained at some 
value h^min > 0. For k> 1, we have u(l) > ri (not shown). 
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The asymptotic expressions in eq. [12] correspond to two different physical regimes (see also the regimes discussed 
in 3): In the range where ~ Ti, the GNET is the NET ti divided by the probability pi = ^ ''^'^ to find the 

ligand in state 1, which indicates a mean- field situation where switching and absorption proceed independently. In 
this range the switching dynamics can be approximated by an effective non-switching diffusion process with diffusion 
constant — Di/pi. In the case of ^ h,li ^ 1, and ui(l) are inversely proportional to the surface of 

the target, similar to the reaction-controlled NET to a partially absorbing hole [IjUJ. This is very different from the 
mean- field situation and indicates the appearance of strong correlations. 
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FIG. 5: Simulation results for the sojourn time mi(1) and the exit time as a function of h for different I2 (marked by 
various symbols). The results are scaled by the NET ri. The continuous line in the left panel is the asymptotic obtained 
from eq. [T2](see also text). The results for u{l) in the right panel are obtained for k, = 0.1. 



To conclude this study of the mean time to activate a target by a switching Brownian ligand, we found the new 
formula eq . [T2l extending the NET to the case of switching. For D2 > Di, switching can significantly decrease the exit 
time compared to non-switching. For small values of li and I2, the mean exit time is very sensitive to changes in these 
parameters (see fig. [S]and fig. [3]), and this property can be exploited to modulate cellular signaling. For example, 
if switching is due to a chemical reaction, the value of li is changed via fci2 by modulating the concentration of the 
reaction partner. 

We finish by giving two examples in signal transduction where switching might be important. The first is related 
to the search time for a promoter DNA-site by a transcription factor (TF), which alternates between a 3-dimensional 
diffusion in the nucleus and one-dimensional diffusion along the DNA ^Q, [23, . We consider a TF diffusing along 
the DNA and switching between two states, triggered by stochastic conformational changes of the TF. In state 1, 
diffusion is slow due to a high affinity for the DNA where the TF carefully scans the DNA base pairs. In state 2, 
diffusion is much faster, but the TF does not accurately scan the DNA. The search time of this can be much reduced 
compared to a TF which carefully checks all the DNA base pairs. The second application illustrates our finding 
that a diffusing ligand can activate a target fast, although it may spend almost all its time in a state where it has 
no affinity for the target. This is relevant if a ligand needs to activate a target in a state where it is also prone to 
some degradation. We found that a ligand can largely avoid degradation and still perform fast target activation by 
switching between two states, such that it spends most of its time in state 2, where it cannot be degraded nor activate 
the target. 
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